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Aims. Giant radio halos are Mpc-scale synchrotron sources detected in a significant fraction of massive and merging galaxy clusters. 
The statistical properties of radio halos can be used to discriminate among various models for the origin of non-thermal particles in 
galaxy clusters. Therefore, theoretical predictions are important as new radio telescopes are about to begin to survey the sky at low 
and high frequencies with unprecedented sensitivity. 
Methods. We carry out Monte Carlo simulations to model the formation and evolution of radio halos in a cosmological framework and 
extend previous calculations based on the hypothesis of turbulent-acceleration. We adopt a phenomenological approach by assuming 
i-C " that radio halos are either generated in turbulent merging clusters, or are purely hadronic sources generated in more relaxed clusters, 

"off-state" halos. 

Results. The models predict that the luminosity function of radio halos at high radio luminosities is dominated by the contribution of 
halos generated in turbulent clusters. The generation of these halos becomes less efficient in less massive systems causing a flattening 
of the luminosity function at lower radio luminosities, as also pointed out in previous studies. However, we find that potentially this 
can be more than compensated for by the intervening contribution of "off-state" halos that dominate at lower radio luminosities. We 
derive the expected number of halos to explore the potential of the EMU+WODAN surveys that will be carried out with ASKAP 
and Aperitif, respectively, in the near future. By restricting to clusters at redshifts < 0.6, we show that the planned EMU+WODAN 
surveys at 1.4 GHz have the potential to detect up to about 200 new radio halos, increasing their number by one order of magnitude. 
A fraction of these sources will be "off-state" halos that should be found at flux level f iA < 10 mJy, presently accessible only to 
deep pointed observations. We also explore the synergy between surveys at different radio frequencies, the Tier 1 LOFAR survey 
£Nj ■ at 150 MHz and the EMU+WODAN surveys at 1.4 GHz. We predict a larger number of radio halos in the LOFAR survey due to 

the high LOFAR sensitivity, but also due to the existence of halos with very steep spectrum that glow up preferentially at lower 
frequencies. These halos are only predicted in the framework of turbulent re-acceleration models and should not have counterparts in 
the EMU+WODAN surveys, thus the combination of the two surveys will test theoretical models. 

Key words. Radiation mechanism: non-thermal - galaxies: clusters: general - radio continuum: general - X-rays: general 

1 . Introduction ters through the acceleration of high-energy particles via shocks 

and turbulence (e.g., Sarazin 2004, Brunetti 201 la, for review). 
Radio halos are diffuse Mpc-scale radio sources observed at the A scen ario proposed to explain the origin of the synchrotron 
H . center of ~ 30% of massive galaxy clusters (e.g., Ferrari et al. emitting electrons in radio halos assumes that seed relativistic 
2008, Ventun 2011, Ferettietal. 2012 for recent reviews). These electrons in the ICM are re-accelerated by the interaction with 
sources emit synchrotron radiation produced by GeV electrons me rger-driven MHD turbulence in merging galaxy clusters (tur- 
diffusing through fiG magnetic fields and provide the most im- bulent re . a cceleration model, e.g., Brunetti et al. 2001 ; Petrosian 
portant evidence of non-thermal components in the intra-cluster 2 001 ). According to this scenario, the formation and evolution of 
medium (ICM). radio halos are tightly connected to the dynamics and evolution 
Clusters hosting radio halos always display evidence of very D f the hosting clusters. The occurrence of radio halos at any red- 
recent or ongoing merger events {e.g., Buote 2001 ; Schuecker et shift depends on the rate of cluster-cluster mergers and on the 
al 2001; Govoni et al. 2004; Venturi et al. 2008; Cassano et al. fraction of the merger energy channelled into MHD turbulence 
2010a). The connection between radio halos and cluster mergers an d re-acceleration of high energy particles. This latter point de- 
suggests that the gravitational process of cluster formation pro- pen ds on microphysics of the ICM and it is difficult to estimate 
vides the energy to generate the non-thermal components in clus- ( se e e.g., Brunetti & Lazarian 201 la). 

Despite model details, one of the most important expecta- 



Send offprint requests to: R.Cassano tions of this scenario is the existence of a population of radio 

* e-mail: rcassano@ira. inaf . it halos with very steep radio spectra which will be mostly vis- 



1 



R. Cassano et al.: Radio Halos in future surveys in the radio continuum 



ible al low radio frequencies (Cassano et al. 2006; Brunetti et 
al. 2008). The existence of this population can be tested with 
the upcoming surveys at low radio frequencies with the Low 
Frequency Array (LOFARfl the Long Wavelength Array (LWA, 
e.g., Ellingson et al. 2009) and the Murchison Widefield Array 
(MWA, e.g., Tingay et al. 2012). Cassano et al. (2010b) have 
shown that the LOFAR Tier-1 "Large Area Survey" at 120-190 
MHz (Rottgering 2010), with an expected rms sensitivity of 
~0.1-0.3 mJy/beam, should detect about 350 giant radio halos 
up to redshift z ~ 0.8, and half of them should have very steep 
radio spectra (with a > 1.9, F{v) oc v~ a ). 

An alternative scenario relies on the production of sec- 
ondary electrons. Inelastic collisions between cosmic ray pro- 
tons and target thermal protons in the ICM produce a popula- 
tion of secondary electrons that can generate diffuse synchrotron 
emission on galaxy cluster-scale {e.g., Dennison 1980; Blasi & 
Colafrancesco 1999). Despite several observations pointed out 
that pure secondary models have difficulties in explaining the 
spectral and morphological properties of several nearby giant ra- 
dio halos {e.g., Brunetti et al. 2008; Donnert et al. 2010; Brown 
& Rudnick 201 1; Jeltema & Profumo 201 1, for additional con- 
straints based on y-ray upper limits), synchrotron emission pro- 
duced by secondary electrons must be present in galaxy clus- 
ters. This comes from the theoretical argument that galaxy clus- 
ters are efficient reservoirs of cosmic ray protons (accelerated 
by structure formation shock waves, injected from radio galax- 
ies, or from supernova driven galactic winds) and consequently 
these protons accumulate with cluster life-time increasing the 
probability to have proton-proton collisions (Volk et al 1996; 
Berezinsky et al 1997; EnBlin et al 1997). 

Recent Fermi-LAT and Cherenkov-telescopes observations 
have placed upper limits on the ratio of non-thermal CRp to 
thermal energy densities at the level of ~ few% in a number of 
nearby clusters (Ackermann et al. 2010; Aleksic et al. 2012). 
Assuming these constraints, it has been suggested that during 
cluster mergers the re-acceleration of secondary electrons by 
compressible MHD turbulence can generate synchrotron radi- 
ation in good agreement with radio halos, while synchrotron 
emission ~ 10 times fainter is generated in dynamically relaxed 
clusters when turbulence is dissipated (Brunetti & Lazarian 
2011b). This theoretical conjecture is consistent with the ob- 
served radio bimodality in galaxy clusters {e.g., Brunetti et al. 
2009; Cassano et al. 2010a), although it predicts a level of emis- 
sion (from pure secondaries) in relaxed clusters that is close the 
upper limits derived for these clusters from present radio obser- 
vations. More recently, Brown et al. (201 1) claimed the detection 
of radio emission from the stacking of SUMMS images of "off- 
state" (non-radio halo) clusters at a level ~ 10 times fainter than 
that of classical radio halos, potentially in line with the above 
theoretical picture. 

In about the next 10 years several revolutionary radio tele- 
scopes will survey the sky with unprecedented sensitivity and 
spatial resolution at very low (LOFAR, LWA, MWA) and 
GHz frequencies (ASKAP; Johnston et al. 2008 and Aperitif; 
Oosterloo, Verheijen, van Cappellen 2010). This gives the oppor- 
tunity to constrain the complex connection between cluster dy- 
namics and diffuse radio emission in galaxy clusters by means of 
statistical studies of adequately large cluster sample. The combi- 
nation of incoming surveys at low (LOFAR) and higher (e.g., the 
EMU survey with ASKAP) frequencies will add considerable 
value in the attempt to discriminate between different physical 
origins of giant radio halos in galaxy clusters and, in general, of 
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non-thermal cluster components. The first LOFAR observation 
at ~ 63 MHz of the giant radio halo in Abell 2256 shows that 
the radio spectrum extracted in the region of the halo steepens at 
lower frequencies. This unexpected result suggests that physical 
scenario more complex than previously thought should be con- 
sidered to explain the formation of the radio halo in this cluster 
(van Weeren et al. 2012). In fact, different populations of rel- 
ativistic electrons may coexist in the volume of the radio halo 
if they originate from different acceleration mechanisms, or if 
electrons are accelerated in a non homogeneous turbulent region 
where the efficiency of particle acceleration change with space 
and time in the emitting volume (van Weeren et al. 2012). 

In this paper, we extend previous statistical modeling of gi- 
ant radio halos by combining a picture based on turbulent re- 
acceleration of relativistic electrons by MHD turbulence in clus- 
ter mergers with the process of continuous injection of secondary 
electrons via p-p collision in the ICM. This provide a novel ap- 
proach to interpret future data from surveys of galaxy clusters. 
This has not the aim to reproduce particular spectral features 
such as those observed in Abell 2256, but its aim is to provide a 
simplified (but viable) description of the general observed prop- 
erties of the population of radio halos in galaxy clusters. 

In Sect. 2 we summarize the main ingredients used in the 
model calculations, derive the occurrence of radio halos in clus- 
ters (Sect. 2.2) and the expected radio halo luminosity functions 
(Sect. 2.3). In Sect. 3 we discuss and model the contribution to 
the diffuse cluster-scale synchrotron emission from secondary 
electrons. In Sect. 4 we describe the EMU and WODAN surveys, 
while in Sect. 5 we derive the expected number of radio halos at 
1 .4 GHz and discuss the potential of the EMU and WODAN sur- 
veys. In Sect. 6 we discuss the potential of combining LOFAR 
and EMU and WODAN surveys in addressing the physics of gi- 
ant radio halos. Our conclusions are given in Sect. 7. 
A ACDM {H = 70 Kms-'Mpc" 1 , Q,„ = 0.3, Q A = 0.7) cos- 
mology is adopted throughout the paper. 



2. Statistical modelling of giant radio halos from 
turbulent re-acceleration 

2. 1 . Main ingredients 

Turbulence generated during cluster mergers may re-accelerate 
relativistic particles and produce diffuse synchrotron emission 
from Mpc regions in galaxy clusters {e.g., Brunetti et al. 2008). 
Diffuse radio emission in the form of giant radio halos should 
be generated in connection with massive mergers and fade away 
as soon as turbulence is dissipated and the emitting electrons 
cool due to radiative losses. Attempts for estimating the statisti- 
cal properties of giant radio halos in the context of this scenario 
have been carried out in the past few years (Cassano & Brunetti 
2005; Cassano et al. 2006; Cassano et al. 2010b). These studies 
allowed for a comparison with the presently observed statisti- 
cal properties of giant radio halos, while at the same time they 
provide predictions to test with future instruments. 

In these calculations we model the properties of the halos 
and their cosmic evolution by means of a Monte Carlo approach, 
taking into account the main processes that play a role in this 
scenario: the rate of cluster-cluster mergers in the Universe and 
their mass ratios, and the fraction of the energy dissipated during 
these mergers that is channelled into MHD turbulence and accel- 
eration of high energy particles. We refer the reader to the papers 
quoted above for details, here for the sake of completeness we 
briefly report the essential steps of those calculations : 
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Fig. 1. Reference spectra of "turbulent" radio halos (solid lines) 
and "off-state" (hadronic) halos (dashed lines) in a massive 
(i.e., M v ~ 2.5 X 10 15 M ; black lines) and less massive 
(i.e., M v ~ 10 15 M ; red lines) cluster. Arrows indicate the posi- 
tion of the steepening frequency, v s , in the two cases. The turbu- 
lent spectra are computed assuming in both cases a merger event 
with a sub-clump of mass AM = 5 x 10 14 M at z = 0.023. 



i) The formation and evolution of dark matter halos of galaxy 
clusters is computed by the extended Press & Schechter ap- 
proach (1974, hearafter PS; Lacey & Cole 1993), which is 
based on the hierarchical theory of cluster formation. Given 
the present-day mass and temperature of the parent clusters, 
the cluster merger history (merger trees) is obtained by using 
Monte Carlo simulations. 

ii) The generation of the turbulence in the ICM is estimated for 
each merger identified in the merger trees. Turbulence is as- 
sumed to be generated (and dissipated) within a timescale of 
the order of the cluster-cluster crossing time in that merger. 
Furthermore, it is assumed that the turbulence is generated 
in the volume swept by the subcluster infalling into the main 
cluster. The injection rate of turbulent fast modes/waves, that 
are used to calculate particle acceleration, is assumed to be a 
fraction, rj,, of the PdV work done by this subcluster. 

Hi) The resulting spectrum of MHD turbulence generated by the 
chain of mergers in any synthetic cluster and its evolution 
with cosmic time is computed by taking into account the 
injection of waves and their damping, due to thermal and 
relativistic particles, in a collisionless plasma. Acceleration 
of particles by this turbulence and their evolution is com- 
puted in connection with the evolution of synthetic clusters 
by solving Fokker- Planck equations and including all the rel- 
evant energy losses of particles. 

iv) Synchrotron losses, particle acceleration and emissivity are 
calculated assuming homogeneous models (see Cassano et 
al. 2006): radio halos are assumed to be homogeneous 
spheres of radius Rh ~ 500/z 5( j kpc, and volume-average 
values of turbulent energy, acceleration rate and magnetic 



field are adopted. We assume a value of the magnetic field 
which scales with the virial mass of clusters, M v as 

<B>=B <M> {—J—) , (1) 

where b > 0, and B<m> is the value of the rms magnetic field 
associated with a cluster with mass < M >^ 1.6 x 10 15 M Q 
(see next Sect, and Sect. 3 in Cassano et al. 2006, for de- 
tails). This scaling is motivated by the results of cosmo- 
logical MHD simulations that found that the magnetic field 
scales with the temperature (and mass) of the simulated clus- 
ters (e.g., Dolag et al. 2002fl 

2.2. Occurrence of Radio Halos with redshift and cluster 
mass 

Stochastic particle acceleration by MHD turbulence is a rather 
inefficient process in the ICM that accelerates electrons at ener- 
gies m e c 2 y max < several GeV, since at higher energies the radi- 
ation losses quench the acceleration process (see e.g., Brunetti 
2011b and references therein). In the case of a homogeneous 
model this implies a steepening of the synchrotron spectrum that 
affects the capability to detect radio halos at frequencies sub- 
stantially larger than the frequency, v s , where the steepening be- 
comes severe. The frequency v s depends on the acceleration ef- 
ficiency , x, and on < B >, as v, cc< B > x 2 K< B > 2 +B 2 cmb ) 2 
(e.g., Cassano et al. 2006, 2010b). Monte Carlo simulations of 
cluster mergers that occur during the hierarchical process of 
cluster formation (Sect. 2.1) allow for evaluating^ from the es- 
timated rate of turbulence-generation and the physical condition 
in the ICM, and consequently to explore the dependence of v s 
on cluster mass, redshift, and merger parameters in a statistical 
sample of synthetic clusters. 

A reference example of spectra of "turbulent" radio halos as 
modeled in the present paper is shown in FigQ] A simplified 
approach to estimate the occurrence of radio halos in surveys 
at different observing frequencies is to assume that only those 
halos with v s > v a can be observable, v being the observing 
frequency. Energy arguments imply that giant radio halos with 
v s > 1 GHz are generated in connection with the most energetic 
merger-events in the Universe. Only these mergers can gener- 
ate enough turbulence on Mpc scales and give the acceleration 
rate that is necessary to maintain the relativistic electrons emit- 
ting at these frequencies (Cassano & Brunetti 2005). Present sur- 
veys carried out at v a ~ 1 GHz detect radio halos only in the 
most massive and merging clusters (e.g., Venturi 2011 and ref- 
erences therein, for a review). The fact that ~ 1/3 of X-ray lu- 
minous galaxy clusters host giant radio halos observed at ~ GHz 
frequencies has been used to constrain the parameters for the 
generation of turbulence in the modeling described in Sect. 2.1 
(rj, as 0.1 - 0.3, Cassano & Brunetti 2005; Cassano et al 2008a). 

In Fig. [2] we plot the fraction of radio halos with v s > 1 GHz 
(black upper line) and the differential contribution to this fraction 
from radio halos with v s in four frequency ranges (see figure 
caption for details). This is obtained by assuming a reference set 
of model parameters, namely < B >- 1.9/iG, b — 1.5, rj t = 0.2 
(see also Cassano et al. 2006). 

At low redshift (z = 0-0.1), the fraction of clusters host- 
ing radio halos with v„ > 2 GHz increases with the cluster mass, 
reaching a maximum value in the case of very massive clusters 



2 Dolag et al. (2002) found a scaling B oc T 2 , which would imply that 
B oc M 133 if the virial scaling M oc T 3/2 is assumed. 
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Fig. 2. Fraction of clusters with radio halos with v s > 1000 MHz (black, upper, solid lines) as a function of the cluster mass in the 
redshift range 0-0.1 (left panel) and 0.4 - 0.5 (right panel). The fractions of clusters with radio halos with v s in different frequency 
ranges are also shown : 1000 < v s < 1400 MHz (red lines), 1400 < v s < 2000 MHz (blue lines), 2000 < v s < 4000 MHz (magenta 
lines), and v > 4000 MHz (green lines). 



(M v > 2 - 3 x 10 15 M ). On the other hand, the fraction of clus- 
ters hosting radio halos with v s in the range 1-2 GHz (red and 
blue lines) reaches a maximum for clusters masses slightly be- 
low ~ 2x 10 15 M Q , and then decreases for more massive systems. 
This behavior is due to the fact that in general, in our model, the 
fraction of clusters with radio halos increases with the cluster 
mass, since more massive clusters are more turbulent, and thus 
are more likely to host a radio halo. The occurrence of halos 
with relatively lower values of v s (v s < 2 GHz) also increase 
with the cluster mass, but lower the value of v s lower the value 
of the mass up to which the occurrence of these halos increases, 
for larger masses their occurrence decrease since those clusters 
preferentially form halos with larger values of v s . 

Consequently we expect that very massive (and hot) clus- 
ters tend to generate giant radio halos with radio spectra flatter 
(higher values of v s ) than those in less massive systems. A ten- 
dency to have halos with flatter spectra in more massive (hot) 
systems have been reported in literature (Feretti et al. 2004; 
Giovannini et al. 2009, Venturi et al. 2012, sub.). If confirmed 
these findings would support our expectations. 

In addition, we find that giant radio halos with higher values 
of v s become rarer with increasing redshift, mainly because of 
the inverse Compton losses that increase with redshift and limit 
the maximum energy of the accelerated electrons. 

2.3. The radio halo luminosity function 

The luminosity functions of radio halos (RHLFs) with v s > Vo 
(i.e., the expected number of halos per comoving volume and 
radio power "observable" at frequency vq) can be estimated by : 



dN H (z) dN H (z) IdPivo) 



dVdP(v ) dMdV dM 



(2) 



10 
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Fig. 3. Radio halo luminosity function at y o = 1400 MHz (black 
lines) for clusters at redshift 0.1 - 0.2. The contributions from 
halos with v s in different frequency ranges are also shown : 
1000 < v s < 1400 MHz (red lines), 1400 < v s < 2000 MHz 
(blue lines), 2000 < v s < 4000 MHz (magenta lines), and 
v > 4000 MHz (green lines). The dashed lines correspond to 
the minimum detectable radio halo power (according to Eq. [9] 
with ^ = 3, see text for details) for EMU + WODAN survey 
and for the NVSS (from left to right). 



where dNnizj/dM dV is the theoretical mass function of radio 
halos with v s > vq, that is obtained by combining Monte Carlo 
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Fig. 4. Radio halo luminosity function at v„=1400 MHz derived 
in the redshift intervals: z=0 - 0. 1 , 0. 1 - 0.2, 0.2 - 0.3, 0.3 - 0.4, 
0.4 - 0.5 and 0.5 - 0.6 (from top to bottom). 



allowing a prompt evaluation of (dP(v )ldM)^ Y ^ in Eq.[3] 
Taking Pi. 4 (v ,M v ) = Pi. 4 (1.4, M,,)(1.4GHz/v ) ff , from Eqg] 
one has: 

i.e., radio halos with synchrotron spectra that steepen at lower 
frequencies will also have smaller monochromatic radio powers 
at the observing frequency v . 

As a relevant example, in Figf3]we report the expected RHLF at 
1.4 GHz (black lines) for z — 0.1 - 0.2 (solid thick lines), where 
we also mark the relative contributions from halos with different 
values of v s (see caption). 

As already discussed in Cassano et al. (2006) and Cassano et 
al. (2010b), the shape of the RHLF flattens at lower radio pow- 
ers because of the expected decrease of the efficiency of particle 
acceleration in the case of less massive clusters. We note that 
halos with v s > 2 GHz (green and magenta lines, Fig.0 do not 
contribute significantly to the RHLF at lower radio luminosities. 
This is because higher-frequency halos are generated in very 
energetic merger events that are extremely rare in smaller sys- 
tems, and have large monochromatic radio luminosities (Eq|5]l 
(red and blue lines, Fig. [3). Finally, we note that the normal- 
ization of the RHLFs decreases with increasing redshift (Fig© 
due to the evolution with z of both the cluster mass function and 
the fraction of galaxy clusters with radio halos (Fig. [2] see also 
Cassano et al. 2006). 



calculations of the fraction of clusters with radio halos and the 
PS mass function of clusters {e.g., Cassano et al. 2006). 

dP(vo)/dM depends on (is proportional to) the unknown 
number of seeds electrons (or that of cosmic ray protons) in the 
ICM. Following Cassano et al.(2006) we estimate dP(vo)/dM 
from the correlation observed for giant radio halos between the 
1.4 GHz radio power, P(1.4), and the mass of the parent clus- 
ters (e.g., Govoni et al. 2001; Cassano et al. 2006). The observed 
correlation allows us to normalize the radio luminosity of gi- 
ant radio halos under the assumption that "all" radio halos with 
v s > 1 GHz follow the correlation. In addition, in the context 
of the modeling described in Sect. 2.1, the value of the deriva- 
tive dP( \ A) I dM depends on the set of parameters (B < m>, b), 
consequently the slope of the observed P(1.4) - M correlation 
(P oc M a " with c?m = 2.9 + 0.4) can be used to constrain model 
parameters selecting an allowed region in the parameter space 
(B < m>, b) (Cassano et al. 2006). In this paper we shall adopt 
a reference set of model parameters, i.e., B<m> = 1.9^G and 
b = 1.5 (and 77, = 0.2), that falls in this region; in this case 
% ~ 3.3. 

To derive the RHLF at frequency v - 1400 MHz, we con- 
sider the contribution of all radio halos with v s > 1 GHz. We 
first obtained the RHLF for halos with v s in the frequency inter- 
val Av !( -, and then combined the different contributions from the 
considered intervals Av J( - : 

dN H (z) _y, dN H (z) ^ , dP(v ) yi 
dVdP(v ) ^dMdV'^-S dM '^1' 

The relation between the monochromatic radio luminosity P(v ) 
of halos with a given v s (with v s > vo) and that of halos with 
v, = 1.4 GHz has been derived by Cassano et al. (2010b) : 

P Vs (v„, M v ) = P Vs (v s , M V )(^J = Pi. 4 (l-4, M v )(^-J , (4) 



3. Secondary electrons 

An additional contribution to the Mpc-scale synchrotron emis- 
sion from galaxy clusters comes from the process of continuous 
injection of secondary electrons via p-p collisions in the ICM. 
A self-consistent modeling of the re-acceleration of primary 
cosmic rays and secondary electrons/positrons by compressible 
MHD turbulence in the ICM has been developed in Brunetti & 
Lazarian (2011b). Implementing their formalism in our cosmo- 
logical Monte Carlo framework is challenging and out of the 
main focus of the present paper. 

Here we adopt a simplified approach based on two separate clus- 
ter radio-populations. We assume that presently observed radio 
halos are mainly generated in merging clusters by particle ac- 
celeration by turbulence. In more relaxed clusters, turbulence 
cannot maintain a population of relativistic electrons emitting 
at the observing frequencies, v a , and in these cases the domi- 
nant contribution to Mpc-scale radio emission is due to the gen- 
eration of secondary particles. As a first approximation we can 
assume that the level of the emission from secondary particles 
in these clusters is stationary, i.e., it does not depend on clus- 
ter dynamics, since the primary protons accumulate in galaxy 
clusters over the cluster lifetime and continuously generate sec- 
ondaries (Blasi 2001)B Under our hypothesis the emission pro- 
duced by secondary particles can be constrained from the lim- 
its derived from deep radio observations of "radio quiet" galaxy 
clusters (e.g., Brunetti et al. 2007, 2009). These limits are about 
one order of magnitude below the radio-X-ray luminosity cor- 
relation for classical giant radio halos (Fig. 0. More recently, 
Brown et al. (2011) have detected diffuse emission from "off- 
state" galaxy clusters by stacking SUMSS images of ~ 100 clus- 
ters. Potentially, this signal can be contributed by secondary par- 

3 We do not consider here possible modifications of the synchrotron 
emission in these clusters due to magnetic field amplification and 
cosmic-ray diffusion (e.g., Kushnir et al. 2009; EnBlin et al. 201 1). 
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tides in less turbulent systems. Motivated by these recent ob- 
servations, here we assume that clusters where turbulence is not 
enough to produce giant radio halos emitting at the observing 
frequency vo (i.e., with v s < vq, see Sect.2) host diffuse radio 
emission powered by pure secondaries with a luminosity that 
is similar to the upper limits in Fig. [5] Following Brown et al. 
(201 1) we refer to these halos as to "off-state" radio halos. 
The massfunction of "off-state" halos is given by: 



dN s H ec (z, v ) dN c ' 
= (l-fm(M,v ))x 



dVdM 



dVdM 



(6) 



where /rh(M, v„) is the fraction of clusters of mass M with ra- 
dio halos due to turbulence re-acceleration (with v s > v a ) and 
dN^/dV dM is the cluster mass function. The luminosity func- 
tion of "off-state" halos is: 

dN s H ec (z,v ) _ dN s H ec {z,v ) dM 
dVdP ~ dVdM * dP ' 



(7) 



We derive dM/dP from the expected relation between the ra- 
dio luminosity of "off-state" halos and the mass (or Lx) of the 
host clusters. In their most simple formulation, secondary mod- 
els predict a correlation between the radio luminosity of halos 
and the cluster X-ray luminosity {e.g., Kushnir, Katz & Waxman 
that is slightly flatter than that of giant radio halos (see 
Brunetti et al. 2009): 



n norm u X A 



B 1 



B 2 + B 2 . 

cmb 



A norm M x 



B 2 



B 2 



cmb 



(8) 



where the second equivalence is obtained by assuming the M v - 
Lx correlation (taken from Cassano et al. 2006). The normaliza- 
tion factors A* orm and A^ orm are derived in order to be consistent 
with the radio upper limits in the plane P\ j - Lx obtained for 
"radio quiet" systems (Brunetti et al. 2007, 2009). Specifically, 
we adopt two different approaches and obtain two scalings (see 
Fig(5]): a) in the first one we use the scaling of the magnetic field 
with the cluster mass in EqQ]with b — 1.5 and B<m> = 1.9//G, 
and adopt P as 4 x 10 23 Watt/Hz for L x = 10 45 erg/sec (blue line 
in Fig(5j); b) in the second one we consider a constant magnetic 
field B = 3//G and a slightly higher normalization, P*5x 10 23 
Watt/Hz for Lx = 10 45 erg/sec (red line in Fig(5]). The latter ap- 
proach maximizes the contribution from secondary electrons. A 
reference example of spectra of "off-state" (hadronic) halos in 
the case a) is show in FigQ] 

Apparently, the adopted scalings for secondaries are not fully 
consistent with non-detections of higher luminosity clusters, 
however some of these clusters are cool-core systems and their 
luminosity is likely to be boosted (up to ~ 2 times) by the emis- 
sion produced in the core (Cassano et al. in, prep). In addition, 
upper limits reported in Fig [5] have been obtained at 610 MHz 
and scaled to 1.4 GHz by adopting a spectral index a = 1.3, 
while the spectrum of secondary halos assumed in the present 
paper is a = 1.0 implying that a fair comparison with theoretical 
models should use shallower limits. 

In Fig. [6] we report the total RHLF, obtained by combining 
the contributions from "turbulent" radio halos and from (purely 
hadronic) "off-state" halos, in the case a) (left panel) and b) 
(right panel). Under our assumptions, "off-state" halos dominate 
the RHLF at lower radio luminosities where the RHLF due to 
turbulent radio halos flattens. 




4 Here we modify Kushnir et al. formalism to account also for the 
case of weak magnetic fields. 
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Fig. 5. Correlation between the radio halo luminosity at 1 .4 GHz 
and the cluster X-ray luminosity. Clusters from the literature 
(filled circle) and clusters from the GMRT RH Survey (Venturi 
et al. 2008; open circles and black arrows) are reported. The red 
crosses are obtained by staking the radio images of clusters from 
the SUMSS survey (Brown et al. 2011). On the same figure we 
also report the two scalings adopted in the present paper for ha- 
los produced by secondary electrons: case a) (blue line) and case 
b) (red line)(see text for details). 



4. The EMU+WODAN Surveys 

4.1. ASKAP: EMU survey 

The Australian SKA Pathfinder (ASKAP) (Johnston et al. 2008) 
is a new radio telescope being built at the Murchison Radio- 
astronomy Observatory in Western Australia. It consists of 36 
12-metre antennas distributed over a region 6 km in diameter. 
ASKAP will have an instantaneous field of view of 30 deg 2 , 
enabling surveys of a scope that cannot be contemplated with 
current-generation telescopes. The ASKAP array configuration 
balances the need for high sensitivity to extended structures with 
the need for high resolution. To achieve this, 30 antennas follow 
a roughly Gaussian distribution with a scale of ~ 700 m, corre- 
sponding to a point spread function of ~ 30" , with a further six 
antennas extending to a maximum baseline of 6 km, correspond- 
ing to a point spread function of ~ 10" . These short spacings of 
ASKAP deliver excellent sensitivity to low-surface brightness 
emission, which is essential for studies of radio emission from 
clusters. 

Science data processing will take place in an automated 
pipeline processor in real time. The on-line imaging uses nei- 
ther natural nor uniform weighting, but instead uses an algo- 
rithm called preconditioning, which, together with multi-scale 
clean, gives a similar sensitivity to uniform weighting at small 
spatial scales, and a similar sensitivity to natural weighting at 
large spatial scales. So near-optimum sensitivity is obtained at all 
scales without needing to reweight the data. The ~ 10/iJy/beam 
rms continuum sensitivity in 12 hours is approximately con- 
stant for beam sizes from 10 to 30 arcsec, then increases to 
~ 20 //Jy/beam for a 90 arcsec beam and ~ 40 //Jy/beam for a 
3 arcmin beam. 
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Fig. 6. Total RHLFs obtained by combining the contributions from "turbulent" radio halos and from (purely hadronic) "off-state" 
halos in the case a) (left panel) and b) (right panel) (see text for details). The RHLF are reported in different redshift interval: 
z=0 - 0.1, 0.1 - 0.2, 0.2 - 0.3, 0.3 - 0.4, 0.4 - 0.5 and 0.5 - 0.6 (from top to bottom). 



One important key project of ASKAP will be EMU, the 
"Evolutionary Map of the Universe" (Norris et al 2011), an all- 
sky continuum survey. The primary goal of EMU is to make 
a deep (10 yuJy/beam rms) radio continuum survey of the en- 
tire Southern Sky, extending as far North as +30 deg. EMU will 
cover roughly the same fraction (75%) of the sky as the bench- 
mark NVSS survey (Condon et al. 1998) , but will be 45 times 
more sensitive, and will have an angular resolution (10 arcsec) 
4.5 times better. Because of the excellent short-spacing uv cov- 
erage of ASKAP, EMU will also have higher sensitivity to ex- 
tended structures such as cluster haloes. 



4.2. APERTIF: WODAN survey 



APERTIF, the new Phased Array Feed (PAF) system that will be 
installed on the Westerbork Synthesis Radio Telescope (WSRT), 
will dramatically increase, at frequencies from 1.0 to 1.7 GHz, 
the instantaneous field of view of the WSRT and its observ- 
ing bandwidth (e.g., Oosterloo et al 2010). Many beams can be 
formed simultaneously for each dish making it possible to image 
an area of about 8 square degree on the sky, which is an increase 
of about a factor 30 compared to the current WSRT. This entire 
field of view will be imaged with 15 arcseconds spatial resolu- 
tion over a bandwidth of 300 MHz with a spectral resolution of 
about 4 km/s. The survey speed of APERTIF, and many of the 
other characteristics, will be very similar to ASKAP. 

The extremely large field of view of APERTIF would enable 
the WODAN (Westerbork Observations of the Deep APERTIF 
Northern-Sky) project (e.g., Rottgering et al. 201 1). This project 
has been proposed with the aim to chart the entire accessible 
northern sky at 1400 MHz down to 10 fiJy rms and about 1000 
deg 2 down to 5 pJy. WODAN will be an important complement 
of the EMU project in the northern sky. 



10 



- I 

. NVSS 


i i 




^ ^ - ' ' 


J 

/ 




/ , 
/ / 


* 


/ / 
/ / 




: / / 

: // 
■ // 




I t i i t I i 



0.05 0.1 0.15 0.2 

z 

Fig. 7. Integrated number of radio halos within a given redshift 
expected in the NVSS follow-up of the XBACS clusters (dashed 
line) compared with the observed values (black points) within 
z=0. 1 and z-0.2 taken from Giovannini et al. (1999). The ex- 
pected numbers of halos are computed by consider / m ,„(z) given 
by Eq. [9] with £i = 3 (short dashed line), and by Eq. Q~T] with 
^2 = 10 (long dashed line). 

5. Number of radio halos in the EMU+WODAN 
survey 

It has been shown that results based on calculations carried out 
according to Sect.2 are consistent with the observed increase of 
the fraction of clusters with radio halos with the cluster mass (or 
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Fig. 8. Minimum flux (left panel) and power (right panel) of detectable radio halos at 1400 MHz obtained by Eq. [9] with = 3 
(solid lines) and by EqfTTIwith £2 = 10 (dashed lines). Calculations are shown for the NVSS (rms=0.45 mJy, Ob = 45 arcsec; black 
upper lines) and for EMU+WODAN (rms=10/Jy, Ob = 25 arcsec; red bottom lines). 



X-ray luminosity; Cassano et al. 2008a) and with the observed 
number of nearby radio halos (Cassano et al. 2006). In particu- 
lar, Cassano et al. (2010b) showed that model expectations also 
produce the flux distribution of giant radio halos observed in 
the redshift ranges 0.044 < z < 0.2 (derived from the NVSS; 
Giovannini et al. 1999) and 0.2 < z < 0.32 (derived from the 
GMRT radio halo survey; Venturi et al. 2007, 2008). 

The EMU+WODAN survey will explore the radio sky with 
a sensitivity ~ 10 times better than present surveys, making it 
possible to test models in a totally unexplored range of radio halo 
luminosities and masses of the host systems. In this Section we 
shall derive the expected number of radio halos at 1400 MHz and 
explore the potential of the upcoming EMU+WODAN survey. 
At this point it is important to estimate the minimum flux of a 
radio halo (integrated over a scale of ~ 1 Mpc) that can be de- 
tected in the survey. We will consider two possible approaches: 
i) a brightness-based criterion and ii) a flux-based criterion. 
;) The criterion based on a threshold in brightness guarantees 
that halos are detected in the images generated by the survey. 
From this threshold we derive the flux (luminosity, z) of ra- 
dio halos that can be detected by assuming a spatial distribu- 
tion of their brightness. The brightness profile of giant radio 
halos is known to smoothly decrease with distance from the 
cluster center {e.g., Murgia et al. 2009) implying that the outer- 
most, low brightness, regions of halos are very difficult to detect. 
However what is important is the capability to detect at least, 
the brightest regions, of the radio halos. Radio halos emit about 
half of their total radio flux within their half radius (Brunetti 
et al. 2007, Fig. 1). Following Cassano et al. (2010b) we esti- 
mate the minimum flux of a halo, / m( „(z), that can be detected 
in the survey by requiring that the mean halo brightness within 
half halo radius (0.5 Oh) is ft times the noise level in the map, 
i.e., 0.5 fmin/Nb ~ ftF ms , where Nb is the number of independent 
beams within 0.50h and F rms is the rms noise per beam. 
This gives: 



f mi „(z) - 1.2X10^; 



/ Fnns V 100 arcsec 2 
llO/iJy/l Oj 



A arcsec 2 / 



where 0h(z) is the angular size of radio halos at a given redshift 
in arcseconds and 0b is the beam angular size of the survey in 
arcseconds. 

ii) A second possible approach to derive /„„„ is to assume that the 
halo is detectable when the integrated flux within 0.5 Oh gives a 
signal to noise ratio £2- This is: 

f min (< 0.5 H ) - ft X Frms, (10) 

From Eq. [10]it follows 

/m ,„(z), 1.43xl0-^ 2 (^VHi^VMl) [mJy] , (11 ) 
UOjuJyA 0b A arcsec/ 

In the following, for a better comparison with previous works 
(Cassano et al. 2010b), we will consider i) as our reference ap- 
proach, and will give some results also based on ii). 

The number of radio halos with flux > f m i n (z) in the redshift 
interval, Az = Z2 — Zi, can be obtained by combining the RHLF 
(dN H (z)/dP(v )dV) and f min (z): 



• Ul, 



dN H (P(v ),z') 
dP(v ) dV 



dP(v„) (12) 



[mJy],(9) 



Estimating ft and £2 is the more critical point in this proce- 
dure. Considering case i), in Cassano et al. (2008a) we analysed 
VLA radio observations in D-array configuration of empty fields 
where we injected fake radio halos in the (u,v) plane of the in- 
terferometric data; the injected radio halos were placed at z=0.1 
assuming a diameter ~1 Mpc. Cassano et al. concluded that ra- 
dio halos become visible in the images as soon as their flux ap- 
proaches that in Eq. [9] with ft ~ 1 — 2. At redshift z > 0.1 the 
number of independent beams that cover the half -radius region 
of the halos decreases, so based on our experience a value £1 ~ 3 
yields a more reliable threshold if we are interested in addressing 
the detection of halos over a wider redshift range, e.g., z=0-0.6. 
Brunetti et al. (2007) and Venturi et al. (2008) injected "fake" ra- 
dio halos in the (u,v) plane of GMRT data to evaluate the flux of 
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detectable radio halos. Based on their findings (Fig. 3 in Brunetti 
et al. 2007) we adopt Eq. [TQwith ^2 - 10. 
As a consistency check we derive the expected number of radio 
halos to be detected in the NVSS up to z < 0.2 and compare them 
with those derived from the cross-correlation of the XBACs sam- 
ple (Ebeling et al. 1996) with the NVSS (from Giovannini et al. 
1999). This is shown in Fig. [7] values of %\ ~3and^2 ~ lOgives 
a predicted number of radio halos consistent with that observed 
in the NVSS0 In the following we will use Eq.|9]with £ 1 = 3 to 
compute the minimum flux of detectable radio halos, and we will 
provide also expectations by considering Eq.Q~T]with £2 = 10 in 
a number of cases. 

Figure. [8] shows / mI „ of giant radio halos as a function 
of redshift (left panel), and the corresponding minimum ra- 
dio luminosity (right panel), obtained according to case i) and 
ii). Calculations are shown for the NVSS (black upper lines) 
and for the EMU+WODAN survey (red lower lines). For the 
EMU+WODAN survey we assume rms — 10/iJy and 0b = 25 
arc sec 0. 

In Fig. [9] (left panel) we show the all-sky number of radio halos 
expected in the EMU+WODAN survey. We consider both giant 
radio halos that originate from turbulent re-acceleration in merg- 
ing clusters and "off-state" halos assuming the optimistic case b) 
in Sect. 3 (see caption for details). We consider the flux limit 
derived according to Eq. [9] with % x = 3 (black lines) and report 
also the total number of halos expected according to Eq.fTTIwith 
£2 = 10 (red lines). 

We note that "off-state" halos are expected to contribute signif- 
icantly (about 30%) to the total number of radio halos that are 
expected in the EMU+WODAN survey. Their contribution is 
larger at lower redshift (z < 0.3). 

We predict that the EMU+WODAN survey will detect up to 100- 
200 radio halos in the redshift range 0-0.6. This will increase the 
number of presently known giant radio halos by about one order 
of magnitude. About 2/3 of these radio halos are expected in 
the redshift range 0.1 - 0.4 (Fig|9j right panel). The number of 
radio halos expected in the EMU+WODAN surveys increases 
by a factor ~ 2 when considering the flux limit given by Eq. Qj] 
This increase is due to an increasing number of "off-state" radio 
halos that become detectable in low redshift (z < 0.2) clusters 
according to this prescription. 

Finally, in Fig[10]we report the all-sky number distribution 
of the radio halos detectable by EMU+WODAN as a function of 
the radio flux. "Off-state" halos contribute potentially at smaller 
radio fluxes, /1.4 < 10 mJy, i.e., at fluxes presently accessible 
only to deep pointed observations. To further highlight the poten- 
tial improvement that will be provided by EMU+WODAN with 
respect to the present statistics of radio halos, in the same figure 
we report the number of radio halos detected in the NVSS survey 
(by inspection of XBAC clusters up to redshift 0.2, Giovannini et 
al 1999) and in the GMRT RH Survey (Venturi et al. 2007, 2008; 
in the redshift range 0.2-0.32). For a sanity/consistency check, 
the number of radio halos in present surveys are also compared 
with model expectations derived for turbulent giant radio halos 
according to Eqsl9f[T2land by taking into account the specifica- 
tion of these surveys. 



5 We note that case i) gives a more fair comparison because 
Giovannini et al. (1999) selected halos and candidate halos from the 
inspection of the NVSS images. 

6 We note that such noise level is of the same order of the confusion 
noise expected in this configuration, however it can be reached after 
subtraction of uv components of the sources detected at more than 100 
cr in the survey images obtained at the highest resolution. 
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Fig. 10. All-sky number distribution of radio halos within z < 
0.6 as a function of the radio flux at 1.4 GHz as expected in 
the EMU+WODAN surveys (blue lines). For comparison, in the 
same figure, the expected distribution of radio halos in the red- 
shift range 0.044 < z < 0.2, compared with that observed in 
the NVSS (black lines and symbols) and in the redshift range 
0.2 < z < 0.32 compared with that observed in the "GMRT 
RH Survey" (red lines and symbols) are also reported. The blue 
dashed line shows the distribution of turbulent generated halos. 



6. Comparison between LOFAR and EMU surveys 

The revolutionary radio telescope LOFAR will carry out surveys 
between 15 MHz and 210 MHz with unprecedented sensitivity 
and spatial resolution (e.g., Rottgering et al. 2006), providing a 
breakthrough in the exploration of the Universe at the low radio 
frequencies. In particular the Tier 1 "Large Area" survey of the 
northern sky is planned to reach sensitivities ~0. 1 mJy/beam in 
the frequency range 120-190 MHz. The combination of observ- 
ing frequency and sensitivity to diffuse emission of Tier 1 make 
this survey the most sensitive survey in the pre-SKA era for the 
exploration of non-thermal radio emission from galaxy clusters 
and large scale structure. 

Based on the hypothesis that giant radio halos originate from tur- 
bulent re-acceleration in merging clusters, Cassano et al. (2010b) 
predict the discovery of about 400 giant radio halos at redshifts 
< 0.6 from the analysis of the Tier 1 survey data. This would 
increase the statistics of these sources by a factor of ~ 20 with 
respect to present day surveys. Remarkably about 1/2 of these 
halos are expected with a synchrotron spectral index between 
250-600 MHz of a > 1 .9, and consequently they would brighten 
only at low frequencies, unaccessible to both present observa- 
tions and future observations with ASKAP and APERTIF. 

The surveys planned with the ASKAP array and with the 
APERTIF onboard of WSRT will provide information comple- 
mentary to that coming from LOFAR. These surveys are planned 
to reach a sensitivity similar to that of the Tier 1 survey in the 
case of extended emission with radio spectral index in the range 
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Fig. 9. Left Panel: expected integral number (all-sky) of radio halos as a function of redshift in EMU+WODAN. Right Panel: 
expected number of RH (all-sky) in redshift intervals in EMU+WODAN. The reference set of model parameters (b = 1.5, B < m> = 
1.9 jt/G, < M >- 1.6 X 1O 15 M and r\, = 0.2) is assumed. In both panels black lines show the total number of halos ("turbulent"+"off- 
state" halos, solid) and of "turbulent" halos (dashed) obtained by considering Eq(9]with i;\ — 3 to derive /„„„(z); while red upper 
lines show the total number of halos ("turbulent"+"off-state" halos) obtained by adopting EqQT|with £2 = 10. Points show the 
integral number of radio halos observed so far. 



a ~ 1 - 1 .2. Their synergy with LOFAR will add considerable 
value to discriminate between different physical origin of giant 
radio halos in galaxy clusters and, in general, to constrain the 
evolution of cluster non-thermal components. 

In this Section we shall focus on the comparison between 
expectations for the statistics of radio halos in LOFAR and 
EMU+WODAN. LOFAR surveys and EMU will overlap by 
about one steradiant, while WODAN will be looking at the same 
sky as LOFAR. 

Here we calculate the expected number of giant radio halos in 
the LOFAR Tier 1 survey following Cassano et al. (2010b). We 
upgrade these calculations (i) by using merger trees from Monte 
Carlo simulations of cluster formation history with increased 
statistics with respect to what previously done (i.e., increased 
number of trials), and (ii) by adopting a tighter grid to sample 
the values of the steepening frequency of radio halos, v s . In ad- 
dition, following Sect. 3 we include the contribution to the halos 
statistics at 120-190 MHz from "off-state" halos in more relaxed 
clusters. We adopt a spectral index a = 1 to estimate the radio 
luminosity of these halos at lower frequencies. 

LOFAR and EMU+WODAN will observe different popula- 
tions of giant radio halos, because LOFAR will explore also the 
population of halos with steepening frequency v s < 1 GHz. As 
a first step we calculate the RHLF of these two population of 
halos. As a relevant example, Fig. QT| shows the RHLF of ha- 
los with v, > 120 MHz (black solid lines) and that of halos with 
v s > 1 GHz (red lines) in two different redshift ranges z=0. 1 -0.2 
(left panel) and z=0.3 - 0.4 (right panel). For a prompt com- 
parison between the two populations, the luminosities of halos 
with v s > 1 GHz, which emit in the frequency band spanned 
by EMU+WODAN, are extrapolated at 120 MHz by assum- 
ing a spectral index a « 1.3. The difference between the two 
RHLF is maximized at higher redshift where the RHLF of high- 
frequency halos shows a dip at radio luminosities « 1 x 10 25 
Watt/Hz (at 120 MHz). It reflects the difficulty of generating 



high-frequency giant halos via turbulent re-acceleration in the 
presence of stronger inverse Compton losses (higher z). Indeed 
those high-frequency halos can be generated at higher redshift 
only by very massive merging events, which however preferen- 
tially produce halos with larger radio luminosities (e.g., Eqs. 4— 
5). On the other hand, we find that radio halos emitting at lower 
frequencies can still be generated efficiently at higher redshifts. 
Consequently the detection at low frequencies of a number of 
radio halos with luminosities ~ 10 25 - 10 26 W/Hz in excess of 
that of radio halos observed at higher frequencies with luminosi- 
ties ~ 5 x 10 23 - 5 x 10 24 W/Hz would confirm our theoretical 
expectations; the sensitivities of the LOFAR Tier 1 survey and 
of the planned EMU survey will be suitable to perform this test 
(Fig. [IB. 

At lower luminosities both the high-frequency and low- 
frequency RHLF are dominated by the contribution of "off- 
state" halos. We note that the number density of these halos at 
lower frequencies is slightly smaller than that at higher frequen- 
cies (red curves are slightly higher than the black curves, Fig.fTTI) 
because at lower frequencies an increasing number of clusters 
can generate giant radio halos via turbulent re-acceleration. In 
other words, some of the clusters hosting low luminosity "off- 
state" ("secondary") halos become "on-state" (turbulent) when 
observed at LOFAR frequencies, and migrate from low to high 
radio luminosities in the LOFAR luminosity function. Fig. QT| 
shows that both LOFAR and EMU+WODAN will be able to de- 
tect "off-state" halos (assuming the optimistic case b) in Sect.3). 

Since LOFAR and WODAN will be looking at the same sky 
with the same sky coverage, here we report a comparison be- 
tween the expected number of radio halos in these two surveys. 
To derive the number of radio halos we follow the procedure in 
Sect. 5, assuming a threshold value £1 ~ 3 for both LOFAR and 
WODAN. In Fig. \l^(left panel) we show the integral number 
of giant radio halos (in the LOFAR sky, i.e., in ~ 20214 square 
degrees) as a function of redshift in the WODAN (red line) and 
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LOFAR (black line) surveys, and in Fig. [T2l( risht panel) we re- 
port the number of radio halos in redshift intervals. The LOFAR 
Tier 1 survey is expected to detect about 500 radio halos, ~ 4 
times more than WODAN. This is because of the better sensi- 
tivity of the LOFAR survey (see also Fig. [TTJ and from the pos- 
sibility to detect low-frequency radio halos {i.e., v s < GHz) at 
the frequencies spanned by LOFAR. The difference between the 
two surveys increases with redshift mainly due to the increas- 
ing number of low-frequency radio halos that are generated in 
cluster mergers at higher redshifts. 

Giant radio halos with very steep spectrum (low-frequency ha- 
los) are only expected in the framework of the turbulent re- 
acceleration models {e.g., Cassano et al. 2006; Brunetti et al. 
2008) thus unveiling a substantial population of these radio 
sources will promptly discriminate among different models pro- 
posed for the origin of diffuse radio emission in galaxy clusters. 

7. Discussion and conclusions 

In this paper we have presented results from Monte Carlo sim- 
ulations to model the formation and evolution of giant radio ha- 
los in the framework of the merger-induced particle acceleration 
scenario (see Sec. 2) and extended previous calculations by in- 
cluding the contribution of secondary electrons generating "off- 
state" halos in more relaxed galaxy clusters (Sect. 3). To com- 
bine these two mechanisms we follow a phenomenological ap- 
proach in which we assume that those clusters where turbulence 
is not sufficient to generate radio halos emitting at the observ- 
ing frequency, v , host "off-state" halos generated only by the 
emission from secondary electrons. We assume that presently 
observed giant radio halos are mainly driven by turbulent re- 
acceleration in merging clusters and constrain the level of the 
"off-state" halos using limits derived for "radio quiet" galaxy 
clusters (Brunetti et al. 2007; Brown et al. 201 1). Under this as- 
sumption "off-state" halos are faint with luminosities typically 
~ 10 times smaller than those of giant radio halos in turbu- 
lent/merging clusters. On the other hand, these "off-state" halos 
are expected to be more numerous than turbulent halos and thus 
they may contribute significantly to the number of radio halos in 
future radio surveys. 

The most important expectation of turbulent re-acceleration 
scenarios is that the synchrotron spectrum of radio halos should 
become gradually steeper above a frequency, v s , that is de- 
termined by the competition between acceleration and energy 
losses and that is connected to the energetics of the merger 
events that generate the halos {e.g., Fujita et al. 2003; Cassano 
& Brunetti 2005). Consequently, in this scenario the population 
of radio halos is expected to be made of a complex mixture of 
sources with different spectra, with massive (and hot) clusters 
that have a tendency to generate halos with spectra flatter than 
those in less massive systems. Contrary to turbulent halos, "off- 
state" halos are expected with power-law spectra with fairly sim- 
ilar slopes. Consequently, surveying the sky at different radio 
frequencies and with appropriate sensitivities allows to disen- 
tangle these two populations. 

In Sects. 2-3 we derive the expected radio halo luminos- 
ity functions (RHLF) at frequency v , that account for radio 
halos originated from turbulent re-acceleration, with steepen- 
ing frequency v s > v,„ and "off-state" radio halos. As a rele- 
vant case we discuss the RHLF at v„ = 1.4 GHz. The RHLF 
is characterized by a flattening at low radio luminosity, P ~ 
few x 10 23 - 10 24 W/Hz, due to the decrease of the efficiency of 
turbulent acceleration in less massive systems, and by an upturn 
at lower radio luminosities due to the contribution of "off-state" 



halos. The flattening/dip is expected to become more prominent 
at higher redshift due to the increase of IC losses that quench the 
acceleration process especially in lower massive systems. 

Future radio surveys have the potential to constrain the for- 
mation and evolution of halos with cosmic time allowing for 
detailed tests of models. Specifically, in Sect. 5 we derive the 
expected number of radio halos in the EMU+WODAN survey. 
The EMU+WODAN surveys will probe the radio sky in the fre- 
quency range 1-2 GHz with a sensitivity 10 times better than 
present surveys. They will allow comparison of model predic- 
tions and observations in a totally unexplored range of radio 
halos luminosities and masses of the hosting clusters. A criti- 
cal point in our paper is to derive a meaningful estimate of the 
sensitivity of these surveys to radio halo emission at different 
redshifts. We explore two possible ways, one based on a bright- 
ness threshold of the radio emission and another one based on 
a threshold in flux density. Threshold values in both cases have 
been estimated from the injection of "fake" radio halos in ex- 
isting surveys (the GMRT RH Survey and the NVSS) under the 
assumption that the brightness distribution of presently known 
halos is representative. A better determination of the sensitiv- 
ity of EMU and WODAN will become available only when 
ASKAP and APERTIF on WSRT will start the commissioning 
phase. Despite the uncertainties on survey sensitivities, the ex- 
pected number of radio halos highlights the potential of the fu- 
ture EMU+WODAN surveys. By using a brightness based cri- 
terion for the detection of halos and assuming the expected sen- 
sitivity of EMU (Norris et al. 2011) and WODAN (Rottgering 
et al. 2011), rms ~ 10yuJy/b, we predict that these surveys will 
potentially discover up to 200 new giant radio halos at redshift 
z < 0.6. Most of these halos are predicted in the redshift range 
z ~ 0.2 - 0.4. This will increase the present number of known ra- 
dio halos by almost a factor 10. The number of halos expected in 
these surveys further increases if a flux based threshold is used 
to estimate the sensitivity of the surveys. In particular, according 
to this method, more radio halos can be discovered at z < 0.2 
with respect to the case of a threshold based on brightness. 
The leap forward of these surveys with respect to present ob- 
servations is due to the fact that the increasing sensitivity will 
make possible the detection of both turbulent radio halos associ- 
ated with less massive systems (or with Lx < 2 x 10 44 erg/s) and 
of "off-state" halos in more relaxed clusters. We derive also the 
flux distribution of expected radio halos, showing that "off-state" 
halos contribute at fluxes /14 < 10 mJy that are presently acces- 
sible only to deep pointed observations; still no clear detection 
of these halos has been obtained so far. 

The most important step in our understanding of the physics 
of radio halos is expected from surveys at lower frequencies and 
from their combination with surveys at higher frequencies. In 
Sect.6 we compare model expectations for the statistics of gi- 
ant radio halos at lower and higher frequencies. We upgraded 
previous calculations used to derive the number of radio ha- 
los at lower frequencies, v a ~ 120 MHz, (from Cassano et 
al. 2010b) i) by using Monte Carlo with improved statistics, 
and ii) by adopting a scenario based on two-populations of ha- 
los (as in Sect. 5, "off-state" and "turbulent" halos). These cal- 
culations allow us to explore the potential of the synergy be- 
tween surveys at lower and higher frequencies with LOFAR 
and EMU/WODAN, respectively. According to the model based 
on turbulent re-acceleration in galaxy clusters different popula- 
tion of giant radio halos should become visible at different fre- 
quencies. Based on our hypothesis of two-populations of halos, 
also the population of "off-state" halos changes with frequen- 
cies because more clusters generate radio halos via turbulent 
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Fig. 11. RHLF of halos with v., > 120 MHz (black solid lines) and that of halos with v s > 1000 MHz (red lines) in two different 
redshift ranges z=0.1 -0.2 (left panel) and z=0.3-0.4 (right panel). As an example, dashed lines (in both panels) show the individual 
contributions of "turbulent" radio halos and "off-state" halos to the LOFAR RHLF. The black and red arrows show the LOFAR and 
EMU+WODAN sensitivities, respectively, at the considered redshifts. 




Fig. 12. Left Panel: integral number of radio halos as a function of redshift in the WODAN (red line) and LOFAR (black line) 
surveys. Right Panel): RH distribution in redshift intervals in the WODAN (red line) and LOFAR (black line) surveys. 



re-acceleration at lower frequencies. We find that the LOFAR 
Her 1 survey will detect about 4 times more halos than the 
WODAN survey, thanks to its better sensitivity and to its lower 
observing frequency that allows the detection of turbulent ra- 
dio halos that become more frequent and luminous at lower fre- 
quencies. The majority of radio halos with very steep spectrum 
(i.e., lower values of v s ) in the LOFAR Tier 1 survey are pre- 
dicted with luminosities P m ~ 10 25 - 10 26 W/Hz at 120 MHz 
and do not have counterparts detectable at higher frequencies. 
As a consequence, the RHLFs at 1 .4 GHz exhibit a dip/flattening 
(at the luminosities P lA ~ 10 25 - 10 26 x (1400/120)" 1 - 3 W/Hz 
~ few x 10 23 - 10 24 W/Hz at 1 .4 GHz). The comparison between 
RHLFs at different frequencies is thus very important as it al- 
lows to promptly unveil the existence of different populations of 



radio halos. For instance, we have shown that comparing for ra- 
dio halos in the LOFAR and in the WODAN surveys will allow a 
prompt test of the existence of radio halos with extremely steep 
spectra. 

We predict that both the LOFAR and EMU+WODAN surveys 
will have the potential to unveil "off-state" radio halos, opening 
the possibility to study diffuse Mpc-scale emission also in more 
relaxed systems. This is currently impossible with present radio 
surveys, and will allow to better understand the effect of clus- 
ter mergers on the evolution of non-thermal radio emission in 
galaxy clusters. 

More generally, LOFAR and EMU+WODAN surveys will al- 
low to readily discriminate between a turbulent re-acceleration 
and a purely "hadronic" origin of radio halos. Indeed, if we as- 
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sume that turbulence does not play a role, and all halos are of 
"hadronic" origin, the luminosity function of "hadronic" halos 
in Fig ITT1 should be boosted up by at least one order of magni- 
tude to explain present number counts, implying a very larger 
number of halos at lower luminosities. Although a quantitative 
comparison between predictions derived by different authours 
is difficult because of the very different hypothesis on (poorly 
constrained) model parameter^ we stress that the most impor- 
tant difference between the two scenarios is in the spectral shape 
of halos and in its consequences. Indeed, re-acceleration models 
lead to the unique predictions of complex spectra, whereas in the 
"hadronic" models the spectra of secondary particles (and thus 
the synchrotron spectra) are rather power-laws which extend in 
principle to very high energies (frequencies). The important con- 
sequences are that in re-acceleration models i) a still hidden pop- 
ulation of radio halos with very steep spectra is predicted to glow 
up at low radio frequencies (e.g., Cassano et al. 2006; Brunetti 
et al. 2008; Cassano et al. 2010b) and ii) the shape of the lu- 
minosity function of giant radio halos depends on the observing 
frequency (FigJTTJ; both predictions will be readily tested by fu- 
ture surveys at different frequencies. 

7.1. Model Simplifications 

Our calculations provide a first step for interpreting future sur- 
veys at the light of models that combine turbulent re-acceleration 
of relativistic particles and the generation of secondary electrons 
in the ICM. The most important simplification in our approach is 
that the evolution of magnetic field in the ICM does not account 
for the possible connection with the level of turbulence, yet the 
magnetic field is simply anchored to cluster mass at any red- 
shift. Present data do not show evidence for a direct connection 
between the magnetic field and the cluster dynamics and turbu- 
lence (Govoni et al. 2010; Bonafede et al. 201 1) thus for the aim 
of the present paper we prefer to keep the models as simple as 
possible. In this paper, we have shown results based on a refer- 
ence set of model parameters (< B >,b and rj,). According to 
Cassano et al. (2006) we expect that the general results given in 
the present paper should be poorly dependent on the adopted pa- 
rameter values. The expected number of radio halos in a given 
survey is expected to change by a factor of about 2 when con- 
sidering different values of parameters (see also Figs. 14 and 15 
in Cassano et al. 2006). Our two-population scenario for halos 
does not consider self-consistently the evolution of protons and 
of their secondaries in the ICM, rather we use a phenomenologi- 
cal approach based on the two "extreme situations" where halos 
originate from turbulence and pure hadronic collisions, respec- 
tively. Brunetti & Lazarian (2011) modeled the acceleration of 
relativistic protons and their secondaries in compressible MHD 
turbulence in a self-consistent way showing that radio halos be- 
come bright when turbulence is generated in the ICM and grad- 
ually evolve into fainter "off-state" halos when turbulence is dis- 
sipated. This complex evolution occurs in a time-scale shorter 
than the life-time of halos and of the hosting clusters, thus our 
approach still provides a valuable way to model the basic statis- 
tical properties of halos (see discussion in Cassano et al. 2010b). 
As a final remark, in this paper we focussed only on the case of 
Mpc-sized radio halos. In reality, smaller halos could be gen- 
erated in dynamically disturbed and relaxed systems by sev- 



7 See, for example, the large discrepancies between Sutter & Ricker 
(201 1) and Zandanel et al. (2012) which attempt to derive expectations 
using the "hadronic" model with different assumptions on the B-mass 
scalings and cosmic ray dynamics. 



eral mechanisms, including turbulent re-acceleration in slosh- 
ing cores and in the region of AGN-driven bubbles (Cassano et 
al. 2008b; Mazzotta & Giacintucci 2008; ZuHone et al. 2012), 
reconnection regions (Lazarian & Brunetti 2011) and hadronic 
collisions (Pfrommer & EnGlin 2004; Keshet & Loeb 2010). All 
these halos constitute an additional population of diffuse syn- 
chrotron source in galaxy cluster to investigate with future radio 
surveys. 

Finally, it's worth mentioning that there are a few radio halos 
found in clusters with X-ray luminosity lower than that of typi- 
cal radio-halo clusters, that are over-luminous in radio (by about 
one order of magnitude) with respect to the radio-X-ray lumi- 
nosity correlation {e.g., Giovannini et al. 2011 and ref. therein). 
These halos have radio luminosities similar to those of classi- 
cal radio halos and are hosted in clusters that are more com- 
mon than the very massive systems hosting classical radio halos. 
Consequently, their observed rarity suggests that they could be 
intrinsically rare. For this reason in our model we do not attempt 
to take into account these sources; future surveys (LOFAR, 
EMU/WODAN) are necessary to get a firm conclusion on the 
occurrence of these sources in clusters. 
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